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The perturbation theory of black holes has been useful recently for providing estimates of gravita- 
tional radiation from black hole collisions. Second order perturbation theory, relatively undeveloped 
until recently, has proved to be important both for providing refined estimates and for indicating 
the range of validity of perturbation theory. Here we review the second order formalism for pertur- 
bations of Schwarzschild spacetimes. The emphasis is on practical methods for carrying out second 
order computations of outgoing radiation. General issues are illustrated throughout with examples 
from "close- limit" results, perturbation calculations in which black holes start from small separation. 



I. INTRODUCTION 



A common occurrence while studying a physical system is that the situation that one is most interested in is too 
complicated to treat in closed form. Yet, for the very same system, it might be possible to study in a closed form a 
situation that has a higher degree of symmetry. An intermediate avenue has usually been to consider the system in 
a regime that is more akin to the one of interest but that departs "slightly" from the situation of high symmetry for 
which one can solve exactly. One is therefore studying "small perturbations" of an exact solution. Considering the 
perturbations as small allows them to be described by a linear theory. The resulting equations usually simplify so 
much that they can be solved in situations without any particular symmetry, provided they depart only slightly from 
an exactly known configuration. 

This approach is applicable in general relativity. The field equations of the theory are a set of non-linear partial 
differential equations that are not solvable in many situations of interest. In fact, state of the art computers might 
not be able to evolve solutions far enough in time for the physics of interest to be captured. Yet, if one assumes 
spherical symmetry, an exact solution has been known since 1916. The exact solution is the Schwarzschild solution 
that describes the geometry of a black hole. As the concept of black hole came to better physical understanding 
in the 1950's, people started studying spacetimes that represented "small departures" from a black hole geometry. 
Regge and Wheeler Q were the first to study these departures, and were able to find a clean formulation for one-half 
of the degrees of freedom of the problem. The formulation for the other half had to wait until the work of Zcrilli 
0. Perturbations of flat spacetime had been considered earlier on ||, and in fact Einstein himself discovered an 
approximate version of the Schwarzschild solution as a perturbation of flat space. 

The motivation of Regge and Wheeler to study perturbations of black holes was to assess the stability of the black 
hole solutions. Were the black hole solutions "an accident" that arose as a consequence of the assumption of spherical 
symmetry only, or did they persist if perturbed? The studies of Regge and Wheeler and Zerilli, and later on of Price 
H and others indeed showed that the black hole solutions are stable under small perturbations. Stability was also 
the motivation of the early second-order perturbations studies due to Tomita and Tajima They were the first to 
workout the second-order perturbation analysis of the Schwarschild solution with the aim of probing the non-linear 
stability of the horizon. 

Since the earliest studies, perturbation techniques have been useful in probing issues of astrophysical significance in 
situations without symmetries, that would otherwise be prohibitively complicated to analyze. Examples of these were 
the analyses of the motion of particles in black hole backgrounds. The formalism was laid out by Zerilli g] and the 
first studies done by Davis et al. 0]. These have subsequently been extended to rotating holes and spinning particles 
(see H and references within). 

The perturbations of the "exterior" Schwarzschild metric could be coupled to perturbations of interior solutions, 
and therefore one could analyze perturbations of stellar objects. This is a subject that lies outside the scope of 
this review. For further references see the work of Cunningham et al. |J and for more recent references see [ jTo| . 
Also outside the scope of this review will be the treatment of perturbations of rotating black holes, first treated by 
Teukolsky [jll| , which very recently has been extended to second order perturbations [Q . 
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The increased activity in numerical relativity applied to black hole collisions |13(] , and the possibility of black hole 
collisions as a source of detectable gravitational waves flifl , has renewed interest in calculations based on perturbation 
theory. Of particular interest have been the "close limit" calculations fL5|] - In this technique, the colliding 

holes start out well inside a single horizon. The initial data, and the spacetime that evolves, outside the horizon, 
are considered to be perturbations of a final stationary hole formed in the collision. Much of this recent work has 
been useful in providing checks of numerical relativity results, providing insight into those results, and in cautiously 
extending some results. The work, so far, has mostly focused on perturbations of Schwarzschild spacetimes, although 
work has very recently started on perturbations of Kerr holes |17|j . 

An inherent feature of linearized perturbation theory, for this or any application, is that there is no "built in" 
indication of how good the perturbation approximation is. For close limit work, linearized perturbation theory 
applies, in principle, only in the limit of zero separation. But in applying this technique, one is not interested so much 
in points of principle as in practical results. One wants to know that a perturbation calculation should agree with 
"perfect" numerical relativity results to, say, within 5%. There are rough indicators of when close limit perturbation 
theory should be applicable. For example, if an apparent horizon surrounds two colliding holes one guesses that close 
limit evolution of the initial data will give tolerable accuracy. In other contexts, there are often similar indicators of 
when perturbation theory applies. But these are rough indicators only, and are typically not quantitative. Short of 
numerical relativity, the only systematic approach to quantifying the errors in linearized theory is the study of higher 
order perturbations. If one computes a physical quantity correctly to second order in a perturbation parameter, then 
the difference between that result, and the result of first order theory alone, is a systematic, quantitative, indication 
of the error in the perturbation theory calculation. 

The need for a good error indicator is paramount in situations where one is "pushing the envelope" of perturbation 
theory, i.e. venturing into large values of the expansion parameter. The case of the "close limit" of black hole collisions 
is one of the first such situations. One is not really interested in the "close limit" of the collision, one is forced into 
it since it is tractable, one is really interested in going as far away as possible from that limit. For that reason a 
significant effort has recently gone into the development of second order perturbations of the Schwarzschild spacetime 
as a practical tool for calculations. We have reported results of second order calculations for head-on collisions of 
equal mass holes Jl^Jl^] , for initially boosted holes [^(| , and for a single slowly spinning hole as it relaxes to its Kerr 
final state fill . Those papers emphasized the results and (especially) the comparison with numerical relativity results. 
In the present paper we supply a more detailed description of how second order perturbation theory is carried out. 
A presentation like this one is important if other researchers are to take advantage of the developed tools in other 
contexts. 

As we mentioned before, there has been some work done in the past on second order perturbations. Tomita and 
Tajima |5| studied second order perturbations in a null formulation, but their interests were in studying the stability 
of the horizon, and therefore their formalism would have to be significantly reworked to address the issue we are 
focusing on, outgoing radiation. Cunningham, et al. || studied the second order perturbations due to the rotation of 
a star, if one viewed the rotation as a perturbation. Because of the narrower range of applicability of these formalisms 
we will not cover them in this review. 

We will, specifically, describe here a scheme for carrying out second order perturbation calculations of the evolution 
of an initial value solution. We will assume that on a t — hypersurface we have a solution of Einstein's initial value 
equations for a 3- metric 7y and an extrinsic curvature K+j, and that these can be expanded to second order in a 
perturbation parameter. We will assume, then, that we have first and second order perturbative initial data. How 
are these data to be evolved forward in time and, in particular, how is outgoing radiation to be computed? 

We will start with a reasonably careful and complete description of first order perturbation calculations. This may 
seem unnecessary in view of the venerable status of first order work. But we will pattern almost all our second order 
equations on their first order equivalents, and it is important that we establish the notation and "style" of the first 
order approach. This will enormously simplify understanding second order calculations since it turns out that almost 
all second order equations are obvious translations of first order equations except for the very important addition of 
"source" terms quadratic in first order perturbations. Patterning the second order equations so closely on the first 
will allow us, when we introduce second order theory in Sec. IV, to focus, not on complexity, but on gauge issues that 
are somewhat subtle, and that have no first order equivalent. These gauge issues are potentially confusing but, as 
will be made clear in Sec. IV, computations of outgoing radiation cannot be carried out without facing these issues. 

The presentation of the first order formalism will be rather general. Prescriptions are given for computations 
starting from any perturbative initial value solution of Einstein's equations, and ending with the outgoing radiation 
amplitudes. We do not give similarly complete and general prescriptions for the second order calculations. The reason 
is the complexity of the source terms that appear in the second order equations. For a second order perturbation 
of a particular multipole index (£, m), and a particular parity (even, odd), there will in general be contributions 
to source terms from first order perturbations of all multipole indices and from both parities. If we were to give 
general source terms, the meaning of the calculations would be obscured by the complexity of the expressions, and 
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the important, but easily overlooked, conceptual issues might be missed. There is a further reason to avoid general 
source term expressions. Such source terms would consist of infinite series of products of first order terms along 
with a set of coefficients in the series. In principle, a notation system could be introduced to represent the Clebsch- 
Gordon-like factors in these coefficients that arise from the projection of a tensor multipole component from the 
product of two tensor multipoles. In practice, such formalism may be unrelated to the practical way in which source 
terms are computed. In examples in which many first order multipoles contribute to each second order multipole, it 
seems plausible that direct numerical methods, rather than formal methods would be used for the computation of the 
source terms. To be concrete, for most second-order applications, the researcher involved will have to construct the 
appropriate source terms. Listing pages and pages of expressions in this paper will not be a viable route to tackle the 
problem, since it is impractical to consider all possible situations. The attempt of this paper is to lay out the general 
formalism of how to deal with higher order perturbations and in particular to address the issues of gauge fixing and 
of meaningful extraction of gravitational waveforms. The lessons developed here are general and useful in all cases. 
The particular expressions will have to be re-worked. 

Clarity, of course, will require examples of explicit expressions. For that reason we will take the following approach. 
Our second order presentation will be restricted to the I = 2, even parity, vacuum perturbations. Generalization to 
other second order multipoles is straightforward, and the equivalent analysis for odd parity is significantly simpler. 
Expressions will not be given for general source terms (source terms for a general system of first order perturbations) 
but will be illustrated with examples from close limit calculations. It turns out that for the initial data sets that have 
been used in close limit work, and until recently as a starting point for most numerical work, the only first order 
perturbations are the £ = 2 even parity perturbations. The Misner initial value solution [2l| ] is a particularly simple 
example of this. In addition to being important, these cases have anomalously simple source terms. This special class 
of examples will be used in Sec. IV to give concrete illustrations of source terms, and more generally of nonlinear 
terms. The reader is cautioned not to mistake these nonlinear expressions as generally valid, . 

The remainder of the paper is organized as follows. Sec. IIA starts with an overview of perturbation theory in 
general including questions of reparameterization and gauge transformations at different orders. The specialization is 
made in Sec. IIB to perturbations of the Schwarzschild spacetime, and the Regge- Wheeler (RW) notation (not to be 
confused with the Regge- Wheeler gauge) is introduced Q . The RW gauge |jj is then introduced, and we emphasize 
the important, but seemingly paradoxical, point that perturbations in the RW gauge can be considered to be gauge 
invariant expressions. First order perturbation theory is given in Sec. Ill, starting in Sec. IIIA, with the presentation 
of the "wave equations," that is, the RW equation and the Zerilli equation g, that are used to evolve perturbations 
in linearized theory. Sec. IIIB discusses the problem of extracting outgoing wave amplitudes from the solution of the 
wave equations, and partially establishes the pattern that will be used for second order work. The second order 
wave equation, a "second order Zerilli equation," is developed in Sec. IVA and is shown, as are all the second order 
equations, to be similar to first order equations except for the inclusion of "source" terms quadratic in first order 
perturbations. The problem of extracting second order outgoing wave amplitude, and the gauge issues involved, are 
taken up in Sec. IVB. Lastly, a summary is given in Sec. V, along with a mention of some related issues. 

A word about our notational conventions is in order. We will use coordinates {t, r, 0, 4>} = {x , a; 1 , a; 2 , a; 3 } always 
to have the meaning of Schwarzschild coordinates, in the limit of small perturbations. Greek indices, a, /3, A, r, /a, v . . . 
will refer to spacetime coordinates, and Latin indices i,j, fc, . . ., will be spatial indices on a constant time surface. 
Throughout the paper we will use a consistent scheme for describing properties of our perturbation functions. A tilde 
(~) denotes that a perturbation quantity is in a general gauge, while "RW" and "AF," as superscripts to the right of 
a symbol, mean that the quantity is defined in the Regge-Wheeler Q or asymptotically flat gauge. Superscripts in 
parenthesis to the left of a quantity will indicate whether the quantity is a first or second order quantity. Superscripts 
in parentheses, to the right of a symbol will (sometimes) be used to denote multipole indices. A numerical subscript 
to the right of a quantity will be used, following the notation of Regge and Wheeler Q], as part of the name of 
various perturbation quantities. To distinguish even and odd parity Regge-Wheeler perturbations, we shall always 
explicitly add a right superscript "odd" to the odd parity perturbations; the omission of the "odd" indicates that the 
perturbation is even. Thus, for example, The quantity ( 2 )^ F ( 3 '°) indicates the axisymmetric octupolar (£ = 3, m = 0) 
part of the second order perturbation of the Regge-Wheeler even parity quantity hi, in an asymptotically flat gauge. 
For compatibility with notation in other papers, and for other purposes, we shall occasionally add further information 
in the form of indices and other attachments to symbols. This makes for notation that is sometimes very cumbersome 
and is often redundant, but we have found this a price well worth paying. In a mathematical description in which first 
and second order quantities, and quantities in different gauges, are used, all with similar symbols, clarity of meaning 
of the symbols is crucial. To reduce the notation we shall usually omit the multipole indices. Of the various indices, 
these seem to us to be the ones that are clearest from context. 

On general notational choices we shall in follow the conventions of Misner, Thorne and Wheeler |2^] and shall use 
the sign convention -H — I — h for the metric, and units in which c = G = 1. A dot over an index shall indicate partial 
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differentiation with respect to t. A subscript following a comma will indicate partial differentiation. For a function of 
a single argument a prime will indicate differentiation with respect to that argument. 



II. PERTURBATION EXPANSIONS IN GENERAL 
A. Basic issues 

Our understanding of perturbation theory in general relativity is based on the idea of g a p(x v ; e) a family of spacetime 
metrics parameterized by some physical quantity e. For our family of spacetimes g a p(x v ] e), we suppose that, at least 
in some restricted region of spacetime, the metric functions can be written as 

g a p{x v ; e) = ^g a p(x") + e Wg a p(x») + e 2 Wg a p(x») (1) 

The metric ' g a js(x u ) is the "background" metric, some known solution of Einstein's equations. The expanded form 
of the metric in (pi) can be substituted in the Einstein equations, and the resulting set of equations can be expanded to 
orders in e. The equations can then in principle be solved order by order, first for W g a p(x"), (first-order perturbation 
theory) then for Wg a p(x v ), and so forth. 

There is an alternative way of viewing first order perturbation theory. One supposes that there is a particular 
spacetime g^p {x v ) that is in some sense close to a known "background" solution (°'g a p of Einstein's equations. One 
writes 

g a p(x»;e) = Wg*0{x v ) + z[glT{x v )- ( °W^)] • (2) 

and substitutes the right hand side of ^ into Einstein's equations, and treats e as a formal parameter for keeping 
track of orders of perturbations. Finally, one sets e to unity. This method turns out to give the correct equations for 
treating g^ (x u ) as a first order perturbation of ^ ' g a p, but it is not based on a systematic approach to perturbation 
theory. For higher order computations, such a systematic approach is important. 

A very relevant example of a family of spacetimes is the Misner spacetimes, representing two equal mass, initially 
stationary throats. The 3-geometry for this initial value solution is given by 

ds 2 = a 2 $ 4 [d^ 2 + drj 2 + sin 2 rj dip 2 ] (3) 

where 

^ y. , =, , * , x • (4) 



-\/cosh(/^ + 2n/io) cos(7y) 



We now introduce new coordinates R,9,(f>, related to fJ*,T),(f> in the same way that spherical polar coordinates are 
related to bispherical coordinates in flat space: 

i2 (R 2 + a 2 ) 2 „ (R 2 -a 2 ) 2 
cosh u = ^ ^ cos 77 = ^5 . (5) 

(R 2 +a 2 f - (2aRcos6) 2 (R 2 - a 2 ) 2 + (2aRsm6) 2 

The metric as presented appears in "isotropic" (conformally flat) form. In order to recover in the close limit the usual 
form of the Schwarzschild solution we transform to a new radial coordinate r through the usual transformation, 

1 



R=-^+V^2Mj , (6) 



where following Misner we have defined M = 4a£i where 



^ — I SmilTl/in 
n— 1 r 



With this new coordinate the 3-geometry takes the form 

dr 

1 ^2M/r 



^Misner = Of f 1 ^TTTZ + ^ d ^ ) > ( 8 ) 
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with 



F=l + 2(1 + ^)- 1 ]T KM(M/Rf +1 P e (co S 9) , (9) 



and with the Kg coefficients given by 



1 ^ (coth n/i )^ , . 

^ o) = ^^g^h^r- ^ 

In the limit fiQ — > 0, all coefficients vanish, so the Misner initial geometry (more properly, that part of it covered 
by r > 2M) approaches the 3-geometry of a constant t slice of the Schwarzschild spacetime. But parameterization 
for perturbation theory of the Misner geometry is nontrivial. The "obvious" choice of parameter fio cannot be 
used, because the deviations from the Schwarzschild geometry are not linear in fi at small [1q. The spacetime 
can be expanded in the parameter <5 = l/ln(/io). The spacetime will then have nonzero perturbations at order 
S , <5 5 , S e , 8 j <5 8 , • • •. If, however, we project out only the quadrupole parts of the deviations from the Schwarzschild 
geometry, the leading terms in the expansion are of order <5 3 , 5 6 , 5 s , S 9 , ■ ■ ■. For consideration of the two lowest orders 
in this expansion it is useful to define e = KzQio) = 0(5 3 ), and to treat the problem as if we were doing perturbation 
to first and second order in e. 



1. Reparameterization 

The expansion in (|l]) is complicated by several issues of choice. One of these is simply the question of the choice of 
parameter j2j|. One could choose a new parameter by e = £[e'], and treat e' as the parameter, so that 

g a p(x";£[e'}) = «V/jCO + ^ ^Va/jCO + ^ ^Ya/jOO • ■ ■ • (H) 

The condition for the spacetime to have a Taylor expansion in both e and e', is that e can be expanded as 

e = Ae' + Be' 2 + ■■■. (12) 

For clarity of explanation, here we will restrict ourselves to transformations with A = 1. (This involves no real loss 
of generality, since it can always be accomplished with a trivial multiplicative rescaling.) With this restriction we 
have g' a B{x v ) = ^ g a p{x v ). The implication is that the first-order prediction for a given numerical value of e (say 
e = 0.1) is identical to the first order prediction for the same value of e' (say e' = 0.1). Since we are free to choose the 
function £, let us suppose that we choose it such that £[0.1] = 10, say from the reparameterization e = e' / (1 + 9.9e'). 
This means that e = 0.1 and e' = 10 correspond to the same physical situation, say the same initial separation between 
two coalescing holes. In the limit as e gets very small, of course, the nature of the transformation in ([l2]) guarantees 
that the first-order perturbation method answers will become insensitive to the choice of parameterization. But that 
fact may be misleading, since in practice one does perturbation calculations for specific physical situations, not as a 
limit. Let us suppose that e = 0.1 represents a physical situation and a parameterization for which first order theory 
gives good results. Then our example of a transformation to e' = 10, leads to inaccurate results despite the fact that 
the underlying physical problem was amenable to a perturbative solution. This illustrates the important point that 
the full nonlinear nature of a choice, like that of parameterization, can have an important effect on the results of first 
order calculations. 



2. Gauge transformations 

Special attention is paid to perturbative coordinate transformations, that is, to transformations that can be written 

as 

x ,f * =x" + e + e 2 (2) f + • • • . (13) 

where C 1 )^, ( 2 )£^, • • • are functions of x^ 1 . The transformation in (|l^) induces a transformation of tensor fields. If 
(suppressing all indices) we let T represent any tensor then the first order transformation is given by 

T' = T-C ( r H T . (14) 
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where L ^ indicates the Lie derivative Q| taken with respect to the vector M£. It is well known that such "gauge" 
transformations play an important role in first order perturbation calculations. Physical answers cannot depend 
on coordinate choices, so the coordinate freedom (gauge freedom) inherent in the gauge transformations must not 
affect physical results. In practice, this can be dealt with by constructing quantities which are invariant under gauge 
transformations, and then using Einstein's equations to compute these quantities. More typically, the problem is 
handled by "gauge fixing," i.e., by imposing specific restrictions on ^g a fi thereby fixing the coordinate system to 
first order. 

For higher order perturbation computations the basic ideas are the same but there are crucial differences in detail. 
First, we remark that the functions ( 2 '£ M in ( |l3| ) are not vector fields, so the transformation induced on tensor 
fields cannot directly be represented by a geometric operation like that in (|l4|). In our second order perturbation 
computations we deal with this by using a two-step process for gauge transformations. 

In the first step the transformation x M = x^ + e is made, and is carried out to (at least) second order. For a 
vector field , for example, the first step is the evaluation of the transformed components by 

V»'(x a ) = (6» + e (1 k» a ) V a {x v ' -e «H • (15) 

The terms on the right can be evaluated exactly, but since we shall end by throwing away everything of order e 3 or 
higher, we need only keep terms of order e 2 . If only these terms are kept the result is 

SV a = V a ' -V a = -e (c (1-)^)" + \<? ~ £2 W r ■ ( 16 ) 

It is of no practical consequence that the e 2 terms do not have a compact geometrical form. It is possible, however, 
to achieve a geometrical form by replacing our simple coordinate transformation x^ — x^ + e W£A< by 

2/ = x" + e Mf* - ±e 2 «£ CT . (17) 

We find that, to second order in e, the induced transformation for any tensor field V is 

SV = -e (C ^V) + \e 2 C (IH (C &6 V) . (18) 

We shall not follow this path here; because of our two step procedure it really makes little difference. What is 
important in the first step is only that all tensor quantities are transformed correctly to second order in e. If the 
first order transformations were not carried out correctly to order e 2 , then our tensor fields would change to that 
order. The spacetime geometry, in particular, would not be isometric to the original spacetime geometry. Though it 
is important that our first order transformation be correct to second order, the nature of the second order coordinate 
change, whether it is that of ([l7|), or the first two terms on the right in (|l^), is irrelevant. The second step will obviate 
any specific choice. 

We call step 1 the "first order gauge transformation" (though it must be carried out correctly to at least second 
order). The purpose of this step will typically be to impose some "first order gauge condition," that is, some conditions 
on the first order metric perturbations ^ ' g a (3 that can be brought about by the proper choice of 

The second step in our procedure is to use a coordinate transformation of the form 

x»' =x" + € 2 {2) C , (19) 
and to transform only to lowest order, i.e., to transform tensors by 

5V = -e 2 (C . (20) 

This transformation does not change the first order parts of tensor fields, and hence leaves intact the first order gauge 
conditions imposed by the transformation in our first step. With the second step we can use the choice of the fields 
( 2 )£ to impose conditions on ^gap- Note that the details of the functions f 2 '^ needed to impose these second order 
gauge restrictions will depend on the second order part of the metric after the first step is performed. The metric, to 
second order is affected by the first step. This means that the details of the second step depend on the details of the 
first step, and this is why we do not need to be specific about the second order part of the coordinate transformation 
used in the first step. 

Our two-step procedure in which we first make the transformation x^ — x^ + e and then x^ = x^ + e 2 ( 2 )£^ 
leads to the following explicit transformations of the metric perturbations. 



G 



'II 
and, 



^ 9uv — ' ^9iiv ^ ^ 9iivip ^ ^ ^9pp ^ \u ' ^ 9pv ^C p ->n (21) 



(2)„' = (2) 



(!)„' (%P_ (!)„' (%P _(!)„' (%J° 

ilpvip s I/pp S !i> J/pj/ S >M 

(0)„ (2)tP_ (0)„ (2)tp _ (0). (2)<rP 



'5p„p (2) e- ( % w (2) ev- (c V (2) ev (22) 

_I (0) o . (%°- (%A _ (0) (1W 

y\i>i<7 s, s, jp .y<rA s 5 p s ^ ■ 

Throughout this paper we shall use the two step procedure for second order gauge transformations. An alternative 
treatment has recently been presented by Bruni et al. |p5fl . That treatment gives a much more geometrical view of 
second order gauge transformations. 



4- Standard hierarchy 

For definiteness we shall now consider specifically the vacuum Einstein equations, which is in fact the case of interest 
for the close-limit method. Einstein's vacuum equations can be written as Q\ T (g a p{ xU )) — 0> where Q\ T represents 
the actions of combining zeroth, first and second derivatives of g a p to form the component G\ T of the Einstein tensor. 
If the expansion in ([l]) is used, the Einstein equations for the family of solutions 

Qxr (g a p(x u ; e)) = Gxr ( {0) g a p{x") + e {1) g a p{x") + e 2 ^ ' g a p{x v ) • • •) - , (23) 

can be expanded in powers of e. The terms to order in e are just those of Q ( ^ g a p{x u )) — 0, and are satisfied 
immediately since, by assumption, ^°'g a p is a solution to Einstein's equations. 
The terms in ( |23|) of first order in e can be written in the form 

eL XT ( Wg afl ) = , (24) 

where L\ T , formally defined by 

(25) 

is a linear operator on ^'g^p, consisting of combinations of differentiation and multiplications by specific coordinate 
functions. The details of L depend on the background solution ^g a p- The equations for ^g a p contained in (S3), 
constitute first order perturbation theory. 

The part of that is proportional to e 2 has terms of two different types. There will be terms linear in ^g a p 
and terms quadratic in ^ g a p- The terms of the first type occur in precisely the same way as do the W g a p terms in 
the first order expression. We can therefore write the second order part of (|||) as 

e 2 L, T ( {2) g a0 )=e 2 ^T XT (Wg a0 ) , (26) 

where ^T\ T is quadratic in W ' g a p. One views (^6|) as a set of linear equations for ^g a p, with the right hand side 
a "source" which is known from the solution of the first order problem. In a similar manner one goes on to find that 
terms in ( |23| ) higher order in e have the form 

e L \ T \9af}) — £ '-i\Ty y 'gap, 'gap, 'g a p,---j ■ {■£<) 

In the n th order source term ^T\ T the combinations of the lower order metric perturbations must occur according 
to obvious rules. For example, the source ^T\ T , for the fifth order perturbations, will in general have contributions 
including terms of fifth power in ^ g a p, terms of second power in ^ g a p multiplied by terms linear in ^ g a p, and so 
forth. (Here "terms in W 'g a p" includes derivatives of these terms.) In principle, one can solve order by order since 



d_ 

de 



Q\t [g<xp{x v ] e)) 
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the source terms for each order are given by the solution known from the lower orders. Furthermore, at each order 
the linear operator being solved is precisely the same L\ T \ all that is changing is the source terms and the initial 
conditions. 

There are alternatives to this "standard hierarchy" of equations for perturbations of increasing order. One could, 
for example, adopt the following alternative iterative scheme. Start with the equation 

G[ (0 W^) + e (1 Wz") + e 2 (2 W^)-") =0. (28) 

As a first step, keep only the terms that are first order in e. This will give a set of linear equations for ^g a /3 identical 
to (g4j). Next, in (p8|), substitute the known solution for ^g a p, omit ^ g a p terms for k > 2 and solve the resulting 
equations to lowest order in ^g a p- The resulting linear equations for ^g a p will not be the same as those in (|26|). 
In the present method we are in effect feeding back all the information about the metric resulting from the first 
order solution, and we are therefore keeping different terms that are higher oder in e. In this alternative method, the 
differential equation for ^g a fi represents propagation of the second order perturbations on a background spacetime 
correct to first order. In the standard hierarchy the second order perturbations, and all perturbations, propagate on 
the zero order background. 

This alternative method would seem intuitively to offer advantages, and indeed is being investigated, but it entails 
a very serious practical difficulty. In the standard hierarchy, the differential operator L embodies the simplicity and 
symmetries of the background solution. In practice this means, for example, that multipole decomposition for angular 
variables can be used. In the alternative method the differential operator for higher order perturbations no longer 
has those simplifications. The higher order equations, as in the standard hierarchy, are linear, but they would need 
to be solved numerically as 2+1 or 3+1 linear hyperbolic systems. 



B. Schwarzschild perturbations 

We specialize now to a spherically symmetric background and introduce coordinates {x°, x , x 2 , x 3 } = {£, r, </>} 
chosen such that in the e -> limit 9, <fi become the usual spherical coordinates. The "background operator" L acting 
on the unknown perturbations embodies the spherical symmetry of the background, and so allows us to eliminate 
angular variables 0, 4>. To take advantage of the symmetry it is necessary to expand the metric perturbations in tensor 
spherical harmonics. This does not require that we make further assumptions about the background geometry, but 
below we shall specialize to the case of a Schwarzschild geometry. To avoid the repetition of very similar lengthy 
expressions, wc give here the description of the multipole decomposition specific to the Schwarzschild background, 
where we follow the notations and conventions |l| of Regge and Wheeler (RW). 

Our background metric ^ g^ v is the standard exterior Schwarzschild solution in Schwarzschild's coordinates 

-(1-2M/V) 

r 2 sin 2 6> 

The 10 metric perturbations can be divided into two sets, called perturbations of odd and even parity, which are not 
mixed by tensorial operators which respect spherical symmetry. The multipole decomposition, is given below for even 
parity perturbations 
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where #o"''"', Hi #2"""', V""', V" , ^' m) and G^.™), are functions of r and t for the even parity parts 
of the perturbations, and (™)/ l ^ ,m ' odd ) j (n)^|X m < odd ) ) (™)ft,^ ,m ^ odd are the odd parity functions. The complete metric 
perturbations are the sum of the even and the odd parity parts. Here we introduce the notation that a tilde ~ over a 
perturbation function indicates that it is in an arbitrary gauge. As an example of metric perturbations we have, from 
(B) and (0), the initial perturbations of the Misner spacetime. For 1 = 1, these are 



R J \ R J 



(32) 



and ^'G — ^'hi — 0. We can choose to set the initial lapse and shift such that ^'Hq = ^'Hi = ^'ho = 0. 

We now fix coordinates to first order by demanding that the perturbations satisfy the RW gauge conditions: For 
even parity this is the condition that the first order even parity functions ^h^ ,m \ ^hf ,m \ Wg^'™) vanish. The 
odd parity RW gauge choice is that the function ( 1 )/ 1 ^' m ) odd vanishes. This specialization is accomplished with the 
first order gauge transformation (pit), with the following notation: 



<V = {1) A Y e m 

= ^A 1 Y i m 

«£2 = WA 2 YJ n ,e+^B odd Y e m , !j> /sm9 

(i)£3 = WA 2 Y e m , 4 ,/sm 2 6- ^ B odA Yp , e / sm.9 



(33) 
(34) 
(35) 
(36) 
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These gauge functions give transformations of the RW metric perturbation functions as follows 

,d WB° dd 
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1 2 <9 WG 
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2 <9r 



Wyli = (1 - 2M/r) 
(D^odd = _j_ (i)^ dd 

With this choice of first order gauge transformation, the result for the first order perturbations is: 

(D /l RW,odd = ^odd + l (i)^ t 

(l^RW.odd = (l)^odd + l r 2 ^-2 (D^° dd ^ r 

and those for the even-parity functions are 

H) K RW = (i)£ +{r _ 2M) ^ (1)^ _ JL (1)^ 

(i) H nw = WH 2 + (2r-3M)( «G, r -A Wfc^ + r(r - 2M) ( ^G, r - \ 



(i) H RW = (DH 1+r 2 (i)g )ir _ (1)^ 



2M 



r(r - 2M) 
2r 



W h + Who, 



r± - 3M) 



2M 



- ^ - M ( W d , r + d,^ + _^_ d,g 



(37) 



(38) 

(39) 

(40) 
(41) 



(42) 
(43) 

(44) 
(45) 
(46) 
(47) 



Here we have dropped the i, m indices in order to simplify the intricate notation. 

The above equations show that we can choose to view the left-hand quantities not as metric perturbations expressed 
in a particular coordinate gauge, but rather (due to the right hand side) as combinations of metric perturbations 
expressed in an arbitrary coordinate gauge. In this sense we can, and we will, view the RW-gauge quantities, such as 
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H RW 

as a compact expression for a general-gauge expression. It should be emphasized that it is trivial, combining 
the above results to construct expressions that are explicitly gauge invariant in terms of the metric perturbations in 
an arbitrary gauge. 

It is instructive to point out the somewhat special aspects of the RW gauge choice and to compare it to other 



approaches. Chandrasekhar |26| chooses to work in a diagonal metric, for axisymmetric perturbations, by making the 
even parity restrictions WiT™ = W/iQ m = ^hf 71 = 0. This choice turns out not to uniquely fix the gauge 
One can, for example, perform a nontrivial transformation ( fIT| ) with W Aq = 0, with ^ A 2 any function of r alone, 
and with ^> A\ given by 

,„ d W A 2 

V>Ai = -r(r-2M) „ 2 ■ (48) 
or 

Such a transformation leaves unchanged the values of ^Hf m , W/iq™ and Wftj m , but changes other perturbations. 
This means that the metric perturbations in the Chandrasekhar gauge are not unique, and thus it is impossible for 
them to be expressed uniquely in terms of arbitrary gauge metric perturbations, as the RW perturbations are in ( |44| ) 

Another approach to perturbation theory is that of Moncrief [|28| who, for analysis of even parity perturbations, uses 

gauge invariant combinations only of the 3-geometry quantities ^H 2 m , ^K e,n , ^G im and WV. This approach 
is quite useful for connecting the computation of radiation to initial value data, since the quantities used are, by 
construction, independent of the (necessarily arbitrary) lapse and shift. The Moncrief approach will be discussed 
again below. Here we confine our attention to the question of whether one could invoke a "Moncrief gauge," a gauge 
in which perturbations of the lapse and shift are set to zero, by choosing (■ 1 >Hf m = ' 1 'if (m = W/iq'™ = 0. In such 
a gauge computations would only involve perturbations of the 3 geometry. It is easy to show, however, that such a 
choice suffers from the same problem as the Chandrasekhar gauge. The gauge conditions do not completely fix the 
gauge. 

The transformation to the RW gauge (or the derivation of the RW arbitrary-gauge expressions) has been discussed 
as a first order problem. We now view this as the first step in the two step process we discussed above. Our second 



step is a second order gauge transformation [see (13)] of the form 

2/ = x» + e 2 (2) ^. (49) 

We choose ( 2 )£ M to impose the RW gauge conditions (e.g., ^G = 0) to second order. The form of ( 2 )£^ needed is 
exactly that of £ M in ( |33| ) through (|4l"|), with the index "1" replaced by "2". Thus for example, second order gauge 
transformations can be written as: 

<V= ^AoYt (50) 
(2)£i = (2) Al ym ( 51 ) 

(2)^2 = m A2 Y i m ) g+^B odd Y e m ^/sme (52) 
(2)^3 = WA 2 Y e m , 4> /sm 2 0- W B odd Y; n , g / sinO . (53) 

To set ^>H[ m = ( 2 )ft,Q m = ^h[ m — we use, for example, ^ A 2 — \ ( 2 'G, and find relations between the second 
order RW perturbations and the second order perturbations in a general gauge, for example 

[2) K RW = (2)K+( r - 2M) ( ^G, r -2r- 2 t 2 ^) . (54) 

It should be mentioned that there appear at second order certain ambiguities in going to the RW gauge if the first 
order perturbations have I — components. The gauge fixing procedure leaves a residual gauge symmetry that has 
to, and can be, dealt with. This issue has been discussed in detail in |2Q], and we will not address it here. 



III. FIRST ORDER PERTURBATIONS 



A. Wave equations 



The first order part of Einstein's equations, the content of (f24|), constitutes linearized perturbation theory. Though 
this theory is long established, we review it here since it forms the foundation of our second order computations. 



11 



Even and odd parity perturbations completely decouple in first order theory and we can analyze them separately. It 
is useful to consider first the relatively simple odd parity problem, first solved in 1957 by Regge and Wheeler In 
the specialized RW gauge the nontrivial odd parity equations are 



d 2 Wfrg w _ d 2 Wfefv _ 2 8 Wfrf w 
dr 2 drdt r dt 



4M 
r 



(1URW 

- (55) 



r(r - 2M) 



-^--a^- + lf - m + 2){r - 2M) -^- = < 56) 

r J dr \ r J dt r z 

Regge and Wheeler [[l] defined, for each multipole, the wavefunction 

godd = r -i (1 _ 2Af/ r )^ W ' odd . (58) 

The field equations above show that this perturbation quantity decouples from ( 1 )/i^ w,odd , and satisfies the wave 
equation 

^-^ + Q odA Vr(r) = 0. (59) 

Here V™ is an ^-dependent function of r that acts as a potential in the wave equation, and r* = r — 2Mln (r/2M — 1) 
is the "tortoise" radial coordinate introduced by Regge and Wheeler. Once (59) is solved and Q odd , and hence / 1 RW ' odd 



are known, one can find /i^ w,odd by solving J57|). In solving that equation, the integration constant is supplied by 
the specification of /i^ W ' odd on the initial t — hypersurface (which is related to the initial extrinsic curvature; see 
below). Once /jR w > odd anc j h^ W odd are known, all gauge invariant odd-parity information is known. 

The odd-parity wave equation (^9|) requires of course the specification of Cauchy data, and these must be supplied 
from the first order perturbations of the 3-metric 7,-j and the extrinsic geometry K%j of an initial value solution. Here 
and throughout we will take our initial surface to be a surface at t — 0. From initial information about the perturbed 

(-\\ flY — '°dd — -odd 

3-metric ( >"fij in some gauge, we can immediately infer the value of ( 'hi at t = and ( >h,2 at t = 0. With 
these, and with p^), we can compute ( 1 )/ 11 RW ' odd at t = and hence can compute <5° dd . The general relationship 
between the extrinsic curvature and the time derivative of the metric is 



Km = — 
3 2N 



d 

9oi\j + 9oj\i - jjjU-j 



(60) 



Here N is the lapse function, given by ^—\/g u . The initial value of the odd-parity shift, ^g°e d -, ^g°$ d or 
equivalently W/iQ dd can be freely specified, and it is convenient to set ( 1 ^g dd initially to zero. With this choice, we 
have 



1 / r d 



Kij= 2\/7^2Mdi 9ij - (61) 

The time derivative of ^h^ 16 follows from the odd parity part of the 96, <jxj>, or 9<ft components of this equation. The 
time derivative of ( 1 '/i° dd follows from the odd parity part of the r9 or rtfi components. From these time derivatives, 
and from @, we find the time derivative of ( 1 )fo 1 RWodd and hence of Q odd . 

Moncricf has taken a distinctly different approach. He works in an arbitrary gauge with only the perturbations of 

, , -odd , . -odd ... -odd 

the 3-gcometry. (Thus for odd parity he works with ( 'hi and ( > h2 , but not ( 'ho .) From those quantities 
he constructs combinations which are invariant with respect to diffeomorphisms on the hypersurfaces. Since the 
quantities, by construction, are automatically independent of shift and lapse choice, they are totally gauge invariant. 
Perturbation quantities which decouple from others must have this property. Moncrief's method leads to the same 
quantity Q odd as that derived by Regge and Wheeler [if Q odd is interpreted in an arbitrary gauge with (|||)]. 

In the case of even parity perturbations, the situation becomes much more complicated. Zerilli has derived a wave 
equation (the "Zerilli equation") for a single decoupled quantity, by working with Fourier transforms, that is, by 
assuming a time dependence e~ luJt . We shall not be able to use Fourier transforms when dealing with the nonlinear 
terms in second order perturbations, so we start by re-deriving the Zerilli wave equation in the time domain. A 
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point to notice is that Zerilli introduces a function R that is equivalent in the time domain to i J Hi dt. To avoid 
introducing an integral we will be working with what amounts to the time derivative of Zerilli's equations. Our time 
domain equations can be compared to the frequency domain equations of Zerilli by replacing d/dt in our expressions 
by — iu>, and by replacing our Wi2 with Zerilli's Rlm, and our ^'x by Zerilli's —i&Klm- 

From one of the vacuum Einstein equations we find that WifJ^ — ^H RW . The remaining equations can be 
broken into a first set that contains only first order derivatives in r and that can be written in the form 
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(62) 
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We also obtain a second set of three equations containing second order derivatives with respect to r. This makes a 
total of six equations for the three unknown functions Wi? RW , WiTfW, and ^K RW , where Wh rw = ^Hf™ = 
^H RW . Compatibility of the system requires then that these equations are not independent. One can show that by 
replacing the first three equations in the second set, one obtains a single compatibility condition (named the "algebraic 
identity" in Ref. M). This can be written as 



£{£+l)M 
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2r 
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Q (1)jjHW 
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g (i) k rw 
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2r 3 d 3 WK KW 



2M 



dt 3 



dt 







(63) 



This equation together with (62), provides a set of four differential equations for ^H RW , ^H RW , and ^K RW . 
One can again show that if (p3) and two of the equations (p2|) are satisfied, then the remaining equation is also 
satisfied. We can now take the t derivative of the equations in fl62|) and use ( p3| ) to eliminate ^H RW , t in the second 
and third of the equations in (62). This reduces the system to two coupled linear partial differential equations for 
(i).ffRW and ( 1 )X RW , with r-dependent coefficients, of first order in r and second order in t. As shown in Zerilli's 
paper this system can be "diagonalized" by the transformation, 



d^K RW /dt = f(r) {1) x + 9{r) (1) R 
«ff RW =h(r) W x + %) {1) R , 



(64) 
(65) 



where 



, A(A + l)r 2 + 3AAfr + 6Af 2 , , 1 
/(r) = r 2 (A.-3Af) ' ^ = 1 



h(r) 



Xr 2 - 'SXMr - 3M 2 



, k(r) = 



A = 



(£-!)(£ + 2) 



(r- 2M)(Xr + 3M) ' " v/ r - 2M ' " 2 
A necessary condition for the Einstein equations to be satisfied is that ^\ and ^'R satisfy 



(66) 



5 (D 



dr* 



X_ (!) 



R , 



dr* 



d 2 

V( r *) + — 
{ dt 2 



(1), 



(67) 



where 



A 2 r 2 [(A + l)r + 3M] + 9M 2 (Ar + M) 
r 3 (Xr + 3M) 2 



(68) 



and, as in (p9|), 



13 



r* =r + 2M ln[r/(2M) - 1] . (69) 
Equations ( |67|) imply that satisfies the "Zerilli equation" 

Q2 (1) Q2 (1) 



dr* 2 dt 2 



V(r*) W X = . (70) 



A necessary condition for perturbations Wflf^ and Wjf- RM '' to satisfy Einstein's equations, is that they satisfy 



the relations in the system of equations (p4|)-(67), for satisfying (JT^) . But in arriving at this system we have 
differentiated Einstein's equations with respect to t. It will be useful in what follows to set down here two of the 
original (not time differentiated) Einstein equations: 

3 B 2 (1) K RW B 2(l) K RW f)(l)K RW 

r(r-2M)i— - ^ + 2(2r-3M) a 



r-2M St 2 v ' dr 2 v y dr 

n (l) ttRW a (l) rrW 

_ 2 (r - 2M)?—!± + 2r , 1 + (£ - l)(t + 2) - 2 . (72) 

or dt 

Once ( [70| ) is solved, the values of all the RW perturbations, can be found. One finds from the first of the 

relations in @, and then from and (|65|) one finds MF RW and «ii: RW . The solution for «ij RW is then 

found from the second of the equations in (|62|). Lastly the "integration constants" in finding Wk rw , WjyKW 

from (D^ RW , «# RW are known from the metric on the t = hypersurface, which fixes the initial values of 
(i)^RW ; (l)^RW 

The Zerilli equation requires Cauchy data, values of ^x and of d ^x/dt, at t = 0, and the Cauchy data must 
originate in an initial value solution of Einstein's equations, i.e., our solution, in some arbitrary gauge, for the 3- 
metric 7y and the extrinsic curvature Kij. From these we immediately get the t = values of W)H 2 , ^'hi, ^'K 
and ^G. To proceed further it is convenient, though not necessary, to choose to view our initial data to be in a 
gauge with (1) h Q = (1) H = = 0. In this gauge (|lj) gives us the time derivatives of (1) # 2 , {1) K, (1) G. 

immediately from Kj. Using this information about the initial 3-geometry perturbations, and the definitions ([44])- 
©, we can find the initial values of Wh rw , ^h rw , Wk rw and of the time derivatives (Vh™, W K RW . 
As an example, from the Misner initial value solution, which has vanishing initial extrinsic curvature, if we choose 
Who = (1) H = {1) Hi = 0, we have 



WiT 2 = W K = 2 



K2 



2M\ 1 /2M\ 3 lA^ 



«G= W hl= CUjj d)G, t = «/M, t = 



(73) 



With these we can solve ( pi]) and ( p5[) for ^X- ^ remains to find ^Xiti which requires 

With the values of the RW perturbations, and derivatives, already found, ^^i7 RW , t follows from ([FT]), and Mif^^t 
from (|72 



The above development of first order perturbation theory gives the time-domain Zerilli formalism in terms of the 
field variable ^x, defined [see ( |64| ) — (p7|)1 by 



(i). 



2M 



Ar + 3M 



r 2 9 (l)^HW 



r-2M 



dt 



(74) 



This development is important here because it will be used as a pattern for the second order formalism. For the first 
order computations themselves, however, it is convenient not actually to use ^x, but rather to use a field variable 
W>tp that is roughly equivalent to the time integral of ^x- 
To introduce we start by using the second relationship in ([32]) to replace ( 1 )/f RW in dTl), with the result 
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(l) x= 



2r{r - 2M) 



£(£+l)(Xr + 3M) 
This leads us to define 

(i)^- 2r(r-2M) 



(i) ijRW _ r ^ _ r-3M (1) ^ RW 



Or 



r - 2M 
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Ar + 3M 
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which, from (f75|) satisfies 



(l)_g-RW 



aWjfBW r _3M (1 



9r 



r-2M 



Ar + 3M 



(i) K RW 



>ijj = *-'x 

It is straightforward to show that ^tp, like ^'X: satisfies the Zerilli equation 

d 2 «V d 2 «V 



dr* 



dt 2 



V(r*) (1 V = 



(75) 



(76) 



(77) 



(78) 



It will be useful for the analysis below of asymptotic behavior to note that in terms of ^ip the equivalents of (|6~j[), 
(f35|), and (|67|) can be written 



\ r J or 



<3r 



1 Wr) + r - 

r / or 



(1) K RW 



(79) 
(80) 
(81) 



The wavefunction ^ip is very closely related to the variable used by Moncrief |p8| , which is defined in terms of 
perturbations in an arbitrary gauge to be 



(1 VMonc 



2r(r - 2M) 
£(£ + l)(Ar + 3M) 



r-2M 



+ 



Ar + 3M 



% + (r-2M) 



<9r r 2 



(82) 



(Note that our expression for d)^ Monc must be multiplied by £(£ + 1) to get the expression in Moncrief 's paper.) It 
is clear that ^'ipMonc reduces to ^ip, when one specializes to the RW gauge and with the use of one of the vacuum 



field equations sets ^'H. 



RW 



(l)^jRW g rea ^. advantage of Moncrief's formulation in (p2[) is that ^^Monc is 



defined entirely in terms of components of the 3-geometry. This greatly simplifies finding Cauchy data from an initial 
value solution. 

The simple relationship of the Zerilli and the Moncrief wavefunctions disappears in the presence of real source 
terms, or of the effective source terms that appear in the second order equations. In our second order development we 
could have followed either the pattern leading to a Zerilli equation for a second order variable analogous to ^x, or 
we could have followed the pattern of Moncrief's derivation p8| to find a Zerilli equation for a variable analogous to 
d^Monc- The former choice has the disadvantage that it is related to the metric perturbations by an additional time 
derivative, but it has the advantage of greater simplicity than the Moncrief approach. The complexity of the second 
order equations makes simplicity the more important consideration, and we have chosen to work with a second order 
equivalent of ^'X- 



B. Radiation 



The solution for first order perturbations does not directly give us the flux of radiated power. To analyze the 
radiation we must go to a coordinate system in which the metric has a manifestly asymptotically flat (AF) form, in 
which the deviations from the Schwarzschild metric fall off with radius as follows. (See Ref. p2|] , Chap. 19.) 
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Sgoo,8g i,Sgu,= 0(r 2 ) Sg Q2 , Sg 03 , Sg 12 , Sg 13 = 0{r 1 ) 

<5#22, <5ff23, ^533 = 0(r) . (83) 

In this coordinate system the information about gravitational radiation is carried by the transverse metric components 
^522, <5ff23, <5<?33, and these perturbations will have the form rx function of t — r* . 

These are requirements on the r dependence of the metric functions for the coordinate system to be AF. If it is to 
be AF to first order then these conditions must be satisfied by the metric perturbations to first order. The RW gauge 
is not AF. There are two rather different ways in which we can extract radiation information from our RW results. 
The first method is to use the gauge invariant expressions for the first order ^'x an d Q odd , and to look at the form 
they take in the radiation zone, for an asymptotically flat gauge. To first order, the AF gauge conditions requires 
that AF perturbations fall of as 

H AF ,H AF ,H AF ~r- 2 h AF ,hf F ,G AF ,K AF -r- 1 (84) 

h AFMd ,h A F ' odd ~ r" 1 h AF > odd ~r (85) 
For first order odd parity perturbations, ([l2|), ( |5^ ) in the radiation zone, then give us, 

godd ^ f _l (l^RW.odd ^ _i r _i (l)/AF,odd ; (86) 

where the dot means differentiation with respect to time, or to retarded time. Of particular interest is the power 
radiated. From the asymptotic form of the perturbations, one finds || 



For even parity, in the radiation zone, using (|74[), the relation equivalent to ( p6|) is: 

(D x « i (r Wk RW - «#f w ) , 
or, using the first order gauge transformation equations, 

and from this asymptotic form, one finds that the radiated power is 

Power ii + 2)[ 1 v 2 (m) 

Power - e i jr^y. WTTr x • (89) 

A very different procedure for extracting radiation is to perform an explicit gauge transformation from the RW 
gauge to a gauge that is AF. Since this will be the basis of the approach we use for the second order problem, we 
illustrate it here. For simplicity of presentation, we limit attention to the even-parity I = 2 case, the case of interest 
for second order calculations. We shall also write the perturbative functions in terms of ^'ip(t, r) instead of ^x(i, r) 
to avoid having to deal with ^K RW , t , and WjJ RW , t , instead of ^K RW , and «7J RW We start by noticing that 
the solution to the Zerilli equation (|T8"|), for outgoing radiation, can be written in the form of an expansion in powers 
of r in which the coefficients are functions of retarded time t — r* . That is, we write d>if)(t, r) in the form 

(l)^( tjr )= «F a (t-r*)+ il) F b (t-r*)/r+ (1) F c (i - r*)/r 2 + 0(l/r 3 ) . (90) 

Substituting this ansatz in ( |78|) one can show that Wfi, (1) F C ,--- are determined once, e.g., ( >F a is fixed. For the 
first terms in the expansion (901), these relations may be conveniently expressed in terms of a suitably defined function 
F(t) as, 

{1) F a (t) = ^ UF»(t) , «F*(t) = , WF,(t) = \ WF(t) - 3 ^ ( ( A A + 1 2 ) ) (1) F'(t) , (91) 

where F'{x) = dF(x)/dx, F"(x) = d 2 F(x)/dx 2 7 etc. In particular, for the quadrupole contribution £ — 2, we have, 
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= I «F"(t -r*) + - {1) F'(t - r *) + \\ «F - M {1) F'] +■■■ . (92) 
3 r r L 

The form of ^tp and hence of W.F(t) will, of course, be determined by the Cauchy data used in the solution of (It). 

For r>Mwe may obtain the asymptotic RW perturbation functions in terms of ^rp, by replacing in ([79]), (8C), 
( jsil ) the asymptotic form of ^ijj. Restricting again for simplicity to £ = 2, we find: 

or r r or 

ttRW ttRW 
#2 — ^0 

From these relations, and from (|9^), we see that the RW perturbations diverge for r — > oo. We can view these 
divergent perturbations as an indication that the RW gauge is not asymptotically flat. (Here we are considering 
the RW quantities as perturbations in a specific gauge, the RW gauge; we usually view them as gauge invariant 
combinations of perturbations) . We can explicitly perform a gauge transformation, analogous to that in (|37|) , to take 
the perturbations from RW gauge to AF gauge: 

(1) rjAF _ (1) „RW , 2M (1) , 2 9 (1)a 

H ° ~ H ° + r(r-2M) Ul + 2 dt 



^H AF = Wi?f w - 



r d W«i r-2M3( 1 'a 



r - 2M dt r dr 



(D,A F = _ r2 ^a + !^( 1 ) ao 

dt r 

«G AF = -2 «a 2 

(Da-af = (DjfRW _ 2 (1)ai 
r 

where we are denoting with ^ao, ^ai, ^a^, the specification of gauge vector the role played by ^Aq, 

^A\, ^^2, in (p?]). To cancel the leading term in ( 1 )_ftT RW , and force the appropriate asymptotic behavior for 
K AF we choose 



rd (1) V> 
2 dr 



= • (95) 



Similarly, the leading divergences in ^Hq W and ^H^ w are cancelled by the choice 



r 



il)a ° = -2^f (96) 
Finally, to avoid a "bad" behavior in M/io and 'hi we choose, 

(1) a 2 = -1 «^ . (97) 
2r 

This procedure may be iterated to the desired degree of accuracy. Since we are only interested in showing that 
there is a choice of gauge functions that carries the RW gauge to an asymptotically flat gauge, the computations are 
simplified if we consider from the start the expression for the RW gauge perturbations in terms of ^F. To put the 
asymptotic form of the full metric in an "asymptotically flat" gauge, we need an expansion of ^'X U P to terms of 
order 1/r 3 . The corresponding asymptotic forms are then, 
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(1) H 



RW 



rd 4 ^F, 2M9 4 Wf, „ 2<9 3 < 1 >F / 



4M 2 d 4 
3 dt 4 
r d 4 Of , 
3 <9i 4 



4M 2 d 4 «F 



RW 



3 <9i 4 

id 3 wf 

3 



HMd 3 d)F „ d 2( -^F, 



7 6 <9i 3 
2M5 4 W.F 



<9i 2 



2 <9 3 W F 



<9i 3 



6 

at 



<9i 3 



<9i 2 



+ 0(l/r 2 ) 



+ 0(l/r 2 ) 



M<9 2 «F, 



1 



2 <9i 2 
0(l/r 2 ) . 



By using these expansions on the right hand side of ( p4| ) we then find for ^ao, and the asymptotic 

expansions, 



(i) 



a Q 



lflWf. 3M9 2 Wf, 2M 2 d 3 «F / 

l m , „ N 3Md^F, „ 3M 2 d 2 WF, 



(t-O 



i 



+ ^A(i-r*) + 0(l/r 4 ) 



(i) 



a i 



4M 3 <9 3 Wf 
~~3 <9i 3 
rd 3 ^F, ts ld 2 ( x )F, . I 

■6-fl?- (t -^-2-^('-^ + r 
1 



laWj? Md 2 ^F, 



+ ^A 1 (t-r*) + 0(l/r 4 ) 



(i) 



i a 2 «f 



a2 



(t-r*) 



6r dt 2 

With these expansions, and (|9q) in 



i a(% i 

(i-r*) + -.4 2 (i-r*) + +0(l/r 4 ) 



1 



3r 2 <9t 
, we finally find 

^M 4 V>F""{t - r*) + ^M 3 W'(t - r*) + 7M 2 «F"(t - r*) 
3 8 4 



+ {1) F{t-r*) + 2A a '{t-r*)+T"{t-r*)] +0(l/r 4 ) 



(DffAF = J_ 



_! M 4 «F""(i - r*) - —M 3 «F"'(f - r*) + -M 2 «F"(i - r*) 
3 v ' 8 v ' 4 1 y 



+ WF(< - r*) - >4o'(* - r*) - -V(< - r*) - T" (t - r*)l + 0(l/r 4 ) 



(99) 



1 
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-=±M 3 V>F"'{t - r*) - -yM 2 V>F"(t - r*) 



+ ( 1 )F(i-r*) + 2^' 1 (i-r*)+F"(t-r*) + 0(l/r 4 ) 
M Q)F"'(t -r*)-- «F'(f - r*) - A' 2 (t - r*) 



12 



0(l/r 2 



M 
12 



V>F"\t -r*)-\ {1) F\t - r*) + A' 2 (t - r*) 



6 



+ 0(l/r 2 ) 



(100) 



(D G af = 1 { i) F , l{ t _ r * } + 2 (1)F , {t _ ^ _ 2_ M{t _ ^ + Q(1/r4) 

(i) k af = 1 (i) F »t t _ r *\ + 2 _ r *) + J_ T M (i) F // t _ r *j + 2 _ r *)l + o(l/r 4 ) 

r r 2r J L 



18 



In all these expressions T is the contribution of order 1/r 4 in an d we notice that Aq,Ai and A2 are not 

determined by the condition of asymptotic flatness, and can be freely specified. 

For computational purposes, it is usef ul to have expressions for the asymptotic forms for Wcto, ^ot\ and 



in terms of ^'ip, instead of ^F as in (100). This inversion is rather cumbersome, because a given term in ^'ip 



contributes to many orders. The leading behavior found in ( |100D is reproduced by, 



2 ot 2 v 

(1) ^ 
2 or at 



m rd (1) ip 1 m d ^ip , . 



but more terms need to b e specified to achieve the same order of accuracy as in (10C). In practice one needs only 
the leading terms in ( |100| ), and having solved for ^ip, one can find ^F, to needed accuracy, from the O(r ) and 
0{r- x ) parts of ^ip. 

We have shown that there are gauge transformations such that the metric functions can be put in the asymptotic 
form, and in which the gravitational wave amplitudes (the O(r ) terms in Wf af and Wg af ) can be found from 
Wip , or, up to an irrelevant, r dependent 'integration constant', in terms of In the next section we explore the 
possibility of obtaining similar results in the second order perturbations. 

IV. SECOND ORDER PERTURBATIONS 
A. Wave Equations 

We now turn to the Einstein equations at second order in the perturbations, and find the first order formalism of 
the preceding section can be modified to handle higher orders. In the second order Einstein equations (^6|), 

Lxr{ i2) g a p) = (2) T AT ( «. 9q/3 ) , (102) 

the left hand side of the second order equations is identical in form to that of the first order equation. The only 
difference is that the perturbations with index n = 1 of first order theory are replaced by those with index n = 2. 
Since the L\ T operator is spherically symmetric for a Schwarzschild background, and linear, L\ T ( ^g a p) can be 
decomposed into multipoles, and the Im multipole of L\ T will involve only that part of ( 2 'g a with the same £m, 
and the same parity. 



The right hand side of (102) can be expanded, as can any functions of {0,<j)), in even and odd multipoles so, for 
example, we can write: 

( 2 )T 00 = ]T WA^Yi (103) 

( 2 )Tno = V ( 2 )«(^. m )£l2L 4- (2)p(l,m)odd 1 01 m 

02 ^ d9 sin6» 



£,m 

and so forth. The multipoles of ^ 2 ^7qo, for example, are found from 



P) A (t,m) = / (2) Too ,yt\» dcos6d(f) . (104) 



It is important to understand that such a multipole decomposition can always be made, but particular im multipole 
of ^T\ T will not be directly related to the same £m multipole of the first order perturbations ^ g a p- The nonlinear 
dependence of ^T\ T on the first order perturbations will mix multipoles and parities. Thus, for example, the I — 2 
second order perturbations will be driven by terms on the right hand side coming from the product of first order 
perturbations with 1 — 2 and those with £ = 4. Even parity second order perturbations will be driven by the products 
of odd parity first order perturbations, etc. 



19 



As an example of this mixing, Cunningham et al. M considered the collapse of a rotating relativistic stellar model, 
to second order in the rate of rotation. The star's first order perturbation is the conserved £ = 1, odd parity, angular 
momentum perturbation. To second order, this rotation drives an even parity I = 2 radiatable perturbation. 

The detailed expressions for the projection of the multipolc components of expressions that are quadratic in mul- 
tipole expansions (as are the ^T\ T ) is tedious, though straightforward B, and will be omitted here. The manner in 
which it is carried out depends on the details of a perturbation problem. In some instances the multipole projections 
might best be carried out numerically. In what follows we shall suppose that the necessary multipole projections have 
been carried out. 

We shall make additional assumptions; we shall consider only the even parity i = 2 second order multipole. 
One justification for this is that it will simplify rather lengthy expressions and greatly simplify their description. 
This will help the presentation to focus on important basic issues about higher order perturbation theory with 
minimal distraction from a minor complication. The generalization from i — 2 to any other multipole is absolutely 
straightforward. We also choose to focus on the quadrupole for a practical reason. Gravitational wave generation seem 
almost always to be dominated by quadrupole radiation, even for sources in which the usual arguments (slow-motion) 
for quadrupole dominance do not apply. The justification for the even parity analysis, is that it is a more difficult 
system to work with. By describing the even parity second order problem we believe we are laying the foundation for 
a reader to make a similar (but simpler) extension from first order to second order odd parity analysis. 

Our next step for second order analysis is to perform a purely second order gauge transformation that sets = 

( 2 )/i x = ( 2 )G = . The second order gauge functions ( 2 )£^ , or equivalently ^Aq, ^A\, ^ A 2 , needed are 
analogous to their first order counterparts, since all changes are taking place purely at the second order. Equations 
(H|) - ( Ejj ) hold true in the second order as well as first; only the indices "(1)" must be changed to "(2)". 

At this point the range of gauge possibilities is wide, and potentially confusing, (i) We could start in some 
arbitrary gauge and transform only the second order perturbations to the RW gauge. That is we could choose to set 

( 2 )/io = ' 2 '/ii = — 0, but not necessarily ^ho — ^hi — ^G — 0. We would then be in a second order, but 
not first order, RW gauge, (ii) We could start from an arbitrary gauge and use a first order transformation (carried 
out to at least second order, of course) to impose the conditions W/io = ^hi = ^G — but not necessarily the 
conditions ( 2 'ho = = ^G = 0. We would then be in a first order RW gauge, but not a second order RW 

gauge, (iii) We could choose to make a first order transformation (carried out at least to second order) to a first order 
RW gauge, followed by a purely second order transformation (which does not affect the first order gauge) to a second 
order RW gauge. We would then have ho = hi = G = to both first and second order, and would be in a first and 
second order RW gauge. 

It is easy to overlook some of the subtleties hidden in the nonlinear interactions of the first and second order 
gauge transformations. One should notice, for example, that the second order gauge functions ^ 2 ^ M used to set 

( 2 ^ho = ^hi — ^G — will depend on whether we have first transformed to the first order RW gauge (be- 
cause such a transformation changes the metric to second order). Since the non- vanishing second order RW parts 

(2) ff RW ; (2)#hw (2)_ ff BW j (2)#BW depend on < 2 ^, we should keep in mind that there are not unique second 
order RW perturbations. The second order RW perturbations, depend on the first order gauge. 



In a second order RW gauge, the second order Einstein equations (102) consist of seven equations linear in the 
second order functions ( 2 )# RW , ( 2 )ff RW , ( 2 )tf 2 RW , and ( 2 )jf RW , but quadratic in the first order functions ^Hf w , 
(1)H$ W , and Wif RW One of these equations is 

(2) jff Rw = WH^ + Sem (105) 

where iSdiff is quadratic in the first order perturbations. Since the first order problem is solved independently, <Sdiff can 
be thought of as a known "source," and we shall use the term "source" to refer below to similar expressions quadratic 
in first order perturbations. 



It is important to understand that the relation (105) exists in a second order RW gauge, whether or not we are 
using a first order RW gauge, but the source term iSdiff will be different (it will have a different numerical value at a 
given coordinate location t, r) depending on the gauge choice that has been made at first order. When it is important 
to emphasize the first order gauge choice that was made in computing source terms like <Sdiff , we w iU use a superscript 
to indicate the first order gauge, so that for example indicates that a first order RW gauge was used, and 
indicates a first order asymptoti cally flat gauge. 

In whatever first order gauge, (|I05|) is the second order equivalent of the first order relationship Wi7 RW = Wff RW , 
and can be used to eliminate ^H^^ in the remaining equations, as was done in the first order case. Note that 
unlike the first order case, we do not now introduce a symbol ( 2 )iJ RW to represent both ( 2 )if RW and ( 2 )ff RW , since 
these second order quantities are not equal. As in the first order case, the Einstein equations consist of two sets of 
three equations for ( 2 'iZ RW , ( 2 'ff RW , and ( 2 )^ RW , one set containing only derivatives of first order in r, and the 
other with derivatives of second order in r. The functions ( 2 '_ff RW , ( 2 'i/ RW , and ( 2 )A" RW appear linearly in these 
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equations, which have the same form, with the same coefficients, as for the system for the corresponding first order 
functions ^H RW , (l)H RW , and ^K RW , but now with "source" terms quadratic in first order perturbations (and 
dependent on the first order gauge choice). The equations of first order in r derivatives are 



Q2 (2) k rw 

drdt 

d 2 ( 2 >ff RW 
drdt 



d ^H™ 
dr 



i a ^H™ 



at 



_ (2) rrRW _ 
„2 n l 



Q2 (2) H RW 



r-3M d^K™ 
r{r - 2M) dt 
AM d^H RW 



S 



K 



-2M dt 2 r(r-2M) 
3M-r d^K RW 3 



r(r - 2M) 
r 



dt 



d ^H RW 



dt 

V>Hf» + S m 

Q (2) K BW 



r - 2M 
-2 



dt 
M 



r{r - 2M) 



r-2M dt 



(106) 



(107) 



(108) 



where Sk and Sh2 are "source terms," quadratic in first order perturbations, and in the r and t derivatives of first 
order perturbations. 

As in the first order case, the second order RW perturbations must satisfy the remaining Einstein equations, three 
equations with second order derivatives in r. These equations are analogous to the corresponding first order equations, 
but the second order equations contain source terms quadratic in the first order perturbations. The procedure that 
was used to simplify the first order system works also in the second order case. For each arrangement of substitutions 
a different source term appears, but one can show that the source terms are equal if the first order perturbations that 
appear in the source terms satisfy the first order perturbation equations. This means that the equations that must 
be solved are three equations of first order in r derivatives, plus an "algebraic identity," of the form 



Q3 (2) K RW Q2 (2) H TM 2j .2 _ 3Mr _ 3M 2 Q (2) R RW g M 



2M 



dt 3 
2r 



dt 2 

SMd^H™ 



r 2 (r-2M) 



dt 



RW 



dt 



S AI = Q 



(109) 



where Sai is a "source" term, quadratic in the first o rder perturbations. In the same ma nner a s in the first order 
case, the next step is to use the "algebraic identity" (109) to eliminate ^H RW :t in (|106|) and (108), and thereby 



to reduce the system to a set of two coupled linear partial differential equations (with "sources") for ( 2 )if RW and 
(2)^rw^ j-j. j s i mm ediate to check that the part linear in ^K RW , t and ( 2 )_ff RW has exactly the same form as the 
corresponding system for W j^ RW and ( 1 )i7 RW in the first order perturbation case, as indeed it must. We may then, 
as in the first order perturbation treatment, introduce a "diagonalization" procedure by defining the functions 
and ( 2 >R, such that 



(2), 



= f(r) {2) X + 9(r)- 



dt 



where f,g,h and k are the same as in 
and ( 2 )if RW is equivalent to 



dt 



^H KW = h(r) ^ X + k{r) 



(110) 



9(2) 



X 



{2) R + Szk 



Using this transformation we may check that the system for ( 2 )i4T RW , t 
d 



dr* ' dr* 

This, again, implies that ^x satisfies an equation of the form 



d 2 

v ^ + dT 2 



(2) 



X + Szr 



dr 



- (2) Y- 
2 X 



dt 2 



(2) 



X~V(r*) + 
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(112) 



where Sz is the "source" term for this second order Ze rilli equation. 

The procedure for specifying the Cauchy data for (112) is patterned on that for first order calculations. From 
our family of initial value solutions, we suppose that we have the 3-metric 7^ and the extrinsic curvature i£y to 
second order as well as first. We assume, for convenience, that the initial value solutions is in a spacetime gauge 
with ( 2 ^ho = ( 2 ">Hi = ^Hq = 0. With (|6l|) we then have the initial values of the perturbation functions 
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( 2 ^G, ^H2, ^hi, and their first time derivatives. We next use the second order equivalent of ( ff4| ) — (p7|) . It 
should be recalled that these equations are the result of a purely second order gauge transformation; the second 
order version of ( ff4| ) — (|47|) is precisely the same as the first order, with "I" superscripts replaced by "2". There 
are no additional terms quadratic in first order perturbations. These relations allow us to compute the initial values 
of ( 2 )i^( RW ), ( 2 )^( RW ) and ( 2 )jj[ RW \ The next step, following the first order pattern, is to use the second order 
Einstein equations 

a (2) ttRW o M fp l {2)jjRW f)(2) K RW\ 

r 9 —^- = — V>Hr + (r - 2M) f ) + S§& (113) 



dt r \ dr dr 



r 3 Q2 (2) K RW Q2 (2) R RW Q (2) R RW 



r(r - 2M) + 2(2r - 3M)- 



r-2Af at 2 v ' 9r 2 v ; dr 

a (2) ff SW Q (2) ttRW 

2( r - 2Af) ^ 2 + 2r ^ + (* - 1)(* + 2) < 2 )^ w - 2 < 2 >tf 2 w + 5™J, (114) 
the second ord er e quivalents of (pfi) and (f72|), to find the initial values of ^ H^- w and ( 2 )_ftT i?w '. With these and the 



definitions in (110) the determination of ^>x an d ^ 2 'X , at t = 0, is complete. 



B. Radiation 

Above we have given a description of how to solve for the second order Zerilli function ^x- It is important to 
understand that this second order perturbation cannot be considered as a second order correction to ^x- Rather, 
we must transform to coordinates which are asymptotically flat (AF) at least to second order in the perturbations. 
In this AF coordinates system, the dominant terms as r — * oo are the transverse K and G terms, which fall off as 
1/r. Most important, in the AF coordinate system the K and G terms give us the intensity of the outgoing radiation. 
The amplitude of the outgoing radiation, correct to first and second order in perturbation theory, therefore, is given 
by K and G correct to first and second order, in a gauge that is AF to first and second order. 

In the problem of finding the radiation to second order there are calculational and conceptual problems that arise 
that are not present in the purely first order problem. In discussing these we start by pointing out how, in principle, 
the second order radiation could be computed in a straightforward way. We could, ab initio, invoke an AF gauge. 
(The specification of a unique AF gauge is not itself well defined, but that is not the issue here.) In this AF gauge 
we could write out the Einstein equations to first and to second order. The source terms ^T a p would have a form 
different from the form we use, based on a ^'T a p computed in a RW gauge. And the various source terms S^f , , 
SjHj, S Af'Sz F would hav e different values than the source terms 5f- w , Sgw S RW > 5 rw 5 rw _ xhe resu lting ^ X 
is a combination of second order RW perturbation functions, and as in the first order case, the relationship to K and 
G in a second order AF gauge can be found by a second order equivalent of a first order procedure. In practice one 
cannot use a first order AF gauge. The enormous complexity of the source term underscores the need for a simple, 
as well as definitive choice of gauge. Our approach is to work in a first order RW gauge and to compute 5 RW , , 

cRW cRW cRW 
°H2 > °AI i °Z ■ 

We will then, in principle compute ^x in a first order RW gauge. In practice, there is a new difficulty not present 
in the first order problem. When the first order RW gauge is used, the source term 6> RW in the second order Zerilli 
equation diverges at large r, and any solution ^x must diverge. This is of course a gauge effect (the effect of 
computing the sources in a first order RW gauge) and not an indication of a physical divergence. But as a practical 
matter we cannot do computations with a divergent quantity. In practice, therefore, we need to modify the problem 
to make it numerically tractable. This is possible since one does not spoil the invariance of ^x under pure second 
order gauge transformations by adding pieces that are quadratic in the first order perturbations. This highlights the 
fact that ^x is f ar from a unique quantity, in fact there is an infinite family of possible candidates for a "second 
order Zerilli function". We therefore add a known term to ^x to cancel the divergent large-r behavior and to define 
a "renormalized" ^'Xn by 

(2) xn= (2) x + r (115) 

which satisfies 
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d 2 < 2 b 



d 2 (2) 



Xn 



dt 2 

d 2 r 



V(r*) ^Xn 



dr* 2 



8 2 T 
W 



(116) 



Here T is a known expression q uadr atic in the (RW gauge) first order perturbations constructed so that the extra 
terms on the right hand side of (jllq) cancel the dominant large-r behavior. Since T is known, numerically solving for 



(2) 



Xn is equivalent to solving for 



2) 



X- 



The explicit forms of source terms depend on the details of the problem and, in particular, on how the various first 
order multipoles contribute. We do not, therefore, give explicit general expressions for the source terms. SW^, <SS^, 



cRW oRW 
°H2 ' °AI > 



S% w , S™- Rather, as an example, we consider the case that contributions to the second order even 



parity quadrupole come exclusively from the first order even parity quadrupole. This specialization, in fact, applies 
to the particular configurations to which second order analysis has already been applied fl£|- ^(J. In this case, the 
rcnormalization is accomplished with 



xl 2) = x (2) 



1 



2r + 3M 



(117) 



and the resulting source iSf'n' , a quadratic expression in the first order perturbation, and some of their space and 
time derivatives, may be written entirely in terms of ^ipi using (|79| ) - (|si|). The explicit expression [] is given in Eq. 
(18),m@. 

— . . (2) 

The process of extracting information about radiation from the solution for \n (t, r) follows the general pattern of 
the first order problem. One starts by writing the solution for ^Xn(t, r) in the form 



(2) 



Xn(t,r)= &>F a (t-r*) + {2) F b (t-r*)/r+ ™F c (t - r*)/r 2 + 0(l/r 3 ) 



(118) 



From the second order Zerilli equation ( |112| ), one can show that ^Fb, ^F c , etc. are determined once ^F a and 
the first order RW perturbations are known. For the special (but important) case that the only relevant first order 
perturbations are the even parity axisymmetric quadrupole perturbations, we have 



(2) 



fi(t) = s| WM' + ^ji 



FfldWFrf 4 ^ d 2 ^Fd 3 ^F 



dt 



dt A 



dt 2 



dt 3 



dt' 3 



dt' 



M 



d 2 (l) Fd 4 (1) F fS^F 



dt 2 dt 4 



dt 3 



(119) 



where is the asymptotic function introduced in (p2|), and we may obtain similar expressions for ^F c (t), and 

higher order coefficients. 

With this form for x( 2 )(2;,r), one then uses the second order system in ( 106 ) — (11C), to solve for the second order 
RW functions. Those functions will, of course, diverge at large r due to now familiar gauge effects. This solution gives 
results, in terms of the asymptotic functions ^F a , ^F b , and ^F c , for the functions ( 2 )if RW , Wh$ w , ^Hf- W , 
( 2 )iJ RW , which are perturbations in a gauge that is RW to both first and second order. To extract information about 
outgoing radiation we must now make two asymptotic gauge transformations. First we must transform from the first 
order RW gauge to an "intermediate" gauge, the result of a first order transformation that makes our coordinates 
system AF to first, but not to second, order. We have already discussed this transformation in Sec. IIIB, in particular 
in equations (|95|) — (101), and found the asymptotic form of the gauge transformation functions ^ag™, W 



Q', 



for each multipole. In Sec. IIIB we were concerned only with the effect of this transformation on first order 



perturbations. Here we are concerned with the effect of the transformations on second order perturbations, and to 



second order, the effect of the transformation will be quadratic in the gauge functions 
will have the form, for example 



and 



( 2 )^ NT = ( 2 > J ff RW + quad, 



(120) 



1 As correctly pointed out by Davies jitj, there is a misprint in reference [^9), the term involving ip 2 should be multiplied by 
(if not it would not be even dimensionally right). For the historical record, the misprint occurred in typing the manuscript, 
i.e. the correct formula was used in the numerical results derived with this formalism, e.g. in reference JL8[. 



23 



where "INT" represents the intermediate gauge, and "quad" stands for some expression quadratic in the first order 
terms, either in products of the ^'a gauge functions, or products of ^'a gauge functions with first order perturba- 
tions. For the first order effects all equations were linear, and each multipole of the Wa!^> m ) gauge functions had an 
effect only on the same multipole of the metric perturbations. In the second order effects, the multipoles mix, and no 
useful general expression can be given for "quad." In our most familiar case, when the only first order perturbations 
of concern are the even parity, axisymmetric, quadrupole perturbations, the explicit gauge transformation is, for 
example, 



(2)^INT = (2) K 



RW 



4 

?7 



-r «a (% RW , t - W>K 



RW 



+9r «a 2 ^K RW -2^ ai 



RW 



-2Mr (1) a 2 - 6r 2 [1) a 2 a 



r 2 (r - 2M) 

+7r 2 «a? + 24r 4 Wa| - 24M 2 Ma§ + 36Mr 2 « ai «a 2 - 18r 3 « ai «a 2 
-A8Mr 3 (1) a 2 2 + 2AMr (1) a 2 -4Mr 2 (1) a 1;r -4Mr 2 (1) a (1) a ut 



(121) 



With such expressions, and with the already established asymptotic forms of the gauge functions given in (99), or 
( |101| ), we now have the relationship between the asymptotic form of an d °f the metric perturbations in the 
intermediate gauge. Equivalently, we have the asymptotic forms of ^Hq NT , ( 2 )i?J NT , ^H™ 1 *, ( 2 )/iq NT , ( 2 )ft,* NT , 
(2) X INT ; (2) G iNT ; in terms of (2)p a (t _ r ^ an d functions related to it. 

The second step of our process is to perform a purely second order transformation to take us from the intermediate 
gauge to a gauge that is AF to both first and second order. We denote the gauge functions that implement this (even 
parity, quadrupolar, axisymmetric) transformation with ^a, that is, our second order gauge transformation uses 



C (2)0 - Wa Q (t,r)Yr , C (2)1 = (2) ai(*,r)17 l 

£(2)2 = (2) u r \^_yra £ (2)3 = (2) n s 1 _g_ ym 

for £, m — 2, 0. These gauge functions, from RW to AF, have the asymptotic form 

Wa = r W a0a (t-r*)+ {2) a ob {t - r*) + ^a 0c (t - r*)/r + • • ■ 
W ai =r &a la (t-r*)+ ^a lb (t-r*)+ (2) a lc (t - r*)/r + • • ■ 
( 2 )a 2 = Wa 2a (t~r*)/r+ ^ a 2b {t - r*)/r 2 + (2) a 2c (< - r*)/r 3 + 



(122) 



(123) 



In the special case that the only first order perturbations that contribute to the second order quadrupole equations 
are the even parity axisymmetric quadrupole perturbations, the coefficient functions are given by 



(2 W(t) 

(2) a la (t) 
{2) a 2a (t) 

(2 W) 

{2) a lb (t) 
{2) a 2b {t) 



WF a (t')dt> - M P>F a (t) - ( (1) F"'(t) 



{2) F a (t>)dt> 
WF a (t")dt" 



(124) 
(125) 
(126) 
(127) 
(128) 
(129) 
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The function is known from the solution to the first order problem, and ^F a is defined in (118) as the 

asymptotic part of the second order Zerilli function. 

The result of the two step gauge transformation is expressions, in terms of ^ F a (t — r*) , and of first order functions, 
for the asymptotic second order metric perturbations in a (first and second order) AF gauge. Among these relations, 
we have ( 2 )G AF and ^>K AF that carry information about the radiation. For the special case of only even parity 
axisymmetric quadrupolc first order perturbations the results is 



d ^G AF 
dt 

d ^k af 

dt 



< 2 >F (t-r*) + 



~ \ {2) F a {t r*) + [ MF»(t r*) MF>"(t - r*) 

2 Fir d 



63 V 5 Si 



WF"(t-r*) V>F"'(t-r*j 



0(r- 2 ) 
0(r- 2 ) 



In terms of the actual results ^'tp and ^ 2 'Xn of first and second order computations, these can be written. 



d (2) G AF _ 1 
dt ~ r 

d ^K AF _ 3 
dt ~ r 



{2) Xn(t,r) + 



2 h_d_ 

5 dt 
2 fn d 



(2 hn{t,r) + '-J'-^( W^f 



7 V 5 dt 



d {1) ^(t,r) 
dt 

d M<ip(t,r) 
di 



0(r- 2 ) 
0(r- 2 ) 



(130) 
(131) 

(132) 
(133) 



All information about gravitational wave energy is carried by G and F in an AF gauge and, in (88), we have seen 
that ( 1 )G AF = ^K AF /3 = ^ 1 'x{^i r )/ r , for £ = 2. We may therefore interpret the expression in brackets in ( 2 )G AF 
as the "second order correction to the gravitational wave amplitude x". In particular, we have that the gravitational 
wave quadrupole power is 



Power 



6?r 
25 



(i) 



(2) 



X + 



tt d 
5 di 



(i) 



ip(t,r) 



d W^t.r) 
di 



(134) 



We may choose to use the expression in (134), or to keep only the terms which are explicitly second order, and 
compute the energy from 



Power = — i Wy • 
25 ' 



2e W; 



(2), 



d ^ip(t,r) 
di 



(135) 



which differs from (134) to third order. For comparison with numerical work |18|j20| | we take the expression in curly 
brackets in (134), aside from overall normalization, to be the gravitational wave amplitude correct to second order. We 
co mpar e radiated energy to the time integral of the power given in (135), although the equally justifiable expression 
in ( 134) turned out to give better agreement with the numerical results used. 



C. Second order techniques for collisions 



The first and simplest application of second order calculations to collisions was the analysis of the collision starting 
with the Misner initial data described in (||) - ( |Io| ) . The family of spacetimes that evolves from these data is described 
by two parameters, but one of them is a trivial overall scaling (say the initial ADM mass). The remaining parameter 
Ho gives a dimensionless measure of initial separation and this parameter (more properly, some function of this 
parameter) is the basis of our ordering of perturbations. [See the discussion following ([To|).] The computations based 
on this scheme have been presented [|l8|, and a comparison given of numerical relativity results, results of first order 
perturbation computations, and the results of perturbation theory to second order. The comparison showed precisely 
the pattern predicted: where the results of first order computations and of second oder computations began to diverge 
(which occurred for parameter /io around 1.8) the results of cither order started to diverge significantly from numerical 
results. This confirmed that second order perturbation would have told us the limiting range of perturbation results 
if we had not had available the results of numerical relativity. 

Subsequent applications of perturbation theory have involved additional complications. Perturbation theory has 
been applied p0| to spacetimes evolving from Bowen-York |nj initial data corresponding to two equal mass holes 
which arc initially moving symmetrically toward each other. In this case we have three parameters. One is an overall 
scaling and can be taken to be the initial ADM mass M of the spacetime. A second parameter is the magnitude 
of the initial momentum P of each hole, and a third is some measure L of the initial distance between the holes. 



25 



Aside from the overall scaling, the family of spacetimes can be characterized by two dimensionless parameters, say 
P/M and L/M. The perturbation analysis reported in |3C]] treated both of these parameters as small. In order to 
apply perturbation methods to such a multiparameter family, it is useful to consider a curve through the parameter 
space. Such a curve gives us a one parameter family of spacetimes, and we can then apply standard methods. The 
curve through the parameter space, however, is not unique. As an example, let us suppose that we are considering 
initial data for P/M = 0.3 and L/M = 0.1. To find perturbation results for this example, we could consider that 
we are on a curve P/M = 3L/M, and L/M is our perturbation parameter. But we could equally well treat the 
spacetime as a point along the curve P/M — 30(L/M ) 2 . The perturbation analysis for close/slow (i.e., small P, small 
L) perturbation theory will depend on which curve in parameter space was chosen. The agreement with numerical 
relativity results will not be equally good for the two choices. 

A different sort of choice was described following (28). It was pointed out that for second order perturbation theory, 
one could "feed back" all information about the first order perturbations; this gives a result that differs only to third 
and higher order from "standard" second order perturbation theory. It was also pointed out that procedure destroys 
the spherical symmetry of the differential operators in the perturbation equations and enormously complicates the 
analysis. 

There is, however, a way of using updating without paying the price of loss of symmetry, and it turns out to give 
an important improvement. We can update only the information about the monopole, thereby preserving spherical 
symmetry of the operators. This turns out to be quite important in the slow/close perturbation analysis of initially 
boosted holes ]3^,Q . The ADM mass of the spacetime turns out to be much more sensitive to the initial momentum, 
than the quadrupole deformation is. A first order change in the monopole mixes with a first order change in the 
quadrupole to produce a second order quadrupole deformation, so our second order radiation computations are 
influenced by the change in the mass. If we use the standard approach, the rapid growth of the ADM mass with 
momentum is the limiting factor that determines the disappointingly small range of applicability of perturbation 
theory. To avoid this limitation we can use a higher order estimate of the ADM mass. We have, in fact, used a 
numerical (rather than perturbative) computation of the ADM mass corresponding to a particular set of parameters 
P and L. This procedure results in an enormous improvement in the range over which perturbation calculations can 
be used for radiation. The reason for this is clear in the following: For values of P large enough so that the ADM 
mass increases by several hundred percent, the radiated energy is still very small. 

This same phenomenon was found in our analysis of the radiation generated as a single spinning Bowen-York ]3l| ] 
hole evolves to its final Kerr state. This problem was analyzed in the limit of slow rotation, so that it could be treated 
as a perturbation of a Schwarzschild hole pq ]. The effect of the spin on the ADM mass was again found to be much 
larger than the effect on radiatable multipoles. Again, a numerical evaluation of the ADM mass was used to extend 
the range of the analysis. 



V. SUMMARY 



Buried, not too deeply we hope, in the equations of the previous section are a few general lessons that justify 
emphasis. We see, for one thing, that in second order work, one relinquishes all the simplicity of linearity that 
characterizes first order perturbation calculations. Gauge transformations, in particular, have a completely different 
character when one works at second order; gauge transformations of first order change second order perturbations. In 
our approach to second order perturbations, these gauge transformations could not be avoided for a rather basic reason: 
One does the mathematics in one gauge, and the physics in others. Actual numerical second order computations have 
to be carried out in a gauge which is both convenient (i.e., minimizes the proliferation of terms) and (unlike the 
"asymptotically flat gauge") is definitive. In our approach that means the use of a system that satisfies the Regge- 
Wheeler gauge conditions to both first and second order. But the Regge- Wheeler gauge is generally not related to two 
other gauges that are of importance to the problem of evolving initial data and finding radiation. First, initial data 
will be found for the problem in some "initial data" gauge, and must be transformed to Regge- Wheeler gauge to give 
us Cauchy data for evolution. Secondly, the result of second order evolution must be transformed, asymptotically, 
to an asymptotically flat gauge in order for information to be extracted about outgoing wave amplitudes. One must 
explicitly perform gauge transformations in the process of making these calculations. And the gauge transformations 
require two steps. One must, in general, perform a first order gauge transformation carried out to second order, and 
next an exclusively second order gauge transformation. 

There are, in principle, other ways to proceed with second order calculations. It is possible, in principle, to construct 
expressions which are formally gauge invariant to both first and second order. The same formal expressions could be 
used for setting the Cauchy data, for evolution, and for interpretation of the results. Such gauge invariant expressions 
would, in effect, have built into them the "gauge invariant" (i.e., unique) character of the Regge- Wheeler gauge. Since 



2G 



the procedure for going from an arbitrary gauge to a (first and second order) Regge- Wheeler gauge is unambiguous 
and local, all second order Regge-Wheeler gauge quantities can be written as combinations of perturbations in a 
general gauge. The development of such a formalism is underway p2|| . By their gauge invariant character such 
expressions would have built into them the gauge transformations, that we perform, in a manner of speaking, as 
an "external" process. The gauge invariant expressions will therefore be significantly more complex in appearance, 
but the important issue is not so much appearance, as suitability for computation. Since these computations (e.g., 
the evaluation of source terms) will be carried out numerically, a gauge invariant expression may turn out be more 
subject to round off error or less, and roundoff error can be a significant problem in expressions combining many high 
derivative terms. 

In closing, we note that second order perturbation theory has turned out to be a great deal more difficult than 
linearized theory, but overcoming these difficulties is motivated by the fact that second order calculations are a great 
deal easier than numerical relativity. What can be (speculatively) said of third and higher order calculations? On the 
one hand we suggest that the step from second order to third and higher order might not be as painful as the step from 
first to second order. The step up to second order required developing the new tools for dealing with nonlinearities. 
With the pattern of those tools established, and with the conceptual issues faced, the next step up adds complexity, 
but, we believe, no new conceptual difficulties. The complexity added by each step up in order is considerable, but 
such complexity is not a crucial obstacle if the work is being done, as it certainly must be, by computers. 

To balance this argument, that one should not be terrified of yet higher order perturbation calculations, it should 
be asked what is to be gained by such higher order results. The overwhelming motivation for second order work 
was the need to establish the range of validity of perturbation analysis. Third order calculations add nothing to 
this, so they are motivated only by the possibility of higher accuracy. But higher accuracy is guaranteed only if the 
perturbation series is convergent. There is no reason to expect that for large values of the expansion parameter the 
series is convergent, and for small values the marginal increase in accuracy would not seem to justify the work of 
going to higher order computations. 
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